
function [wage] = w(det_earn, i,j,k, dt )

wage=exp(det_earn(1) + det_earn(2)*i + det_earn(3)*i*i + det_earn(4)*i*i*i + det_earn(5)*j + det_earn(6)*j*j + k)*dt; 

end
